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CN) ' We compare the rise times of nearby and distant Type la supernovae (SNe la) as 

a test for evolution using 73 high-redshift spectroscopically-confirmed SNe la from the 
J> ■ first 2 years of the 5 year Supernova Legacy Survey (SNLS) and published observations 

of nearby SNe. Because of the "rolling" search nature of the SNLS, our measurement is 

£f) , approximately 6 times more precise than previous studies, allowing for a more sensitive 

test of evolution between nearby and distant SNe. Adopting a simple t 2 early-time 
model (as in previous studies), we find that the rest-frame B rise times for a fiducial 
SN la at high and low redshift are consistent, with values 19.10^q'} 8 (stat) ± 0.2 (syst) 

^' 
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and 19.58l£lg days, respectively; the statistical significance of this difference is only 1.4 
<r. The errors represent the uncertainty in the mean rather than any variation between 
individual SNe. We also compare subsets of our high-redshift data set based on decline 
rate, host galaxy star formation rate, and redshift, finding no substantive evidence for 
any subsample dependence. 

Subject headings: cosmology: observations — supernovae: general 

1. INTRODUCTION 

Type la supernovae (SNe la) have come to play a critical role in attempts to pin down the 
cosmological parameters. Their utility arises because they appear to constitute a class of high- 
quality standardizable candles. Considerable effort has been devoted by many groups to testing 
this assumption. Perhaps the most pernicious concern is that the properties of SNe la have evolved 
between the current epoch and z ~ 0.3 — 1, which characterizes most current "distant" supernova 
samples. One test for evolutionary effects is to compare the rise time, the time from explosion 
until maximum luminosity, of the nearby and distant samples. The early light-curve behavior 
should be influenced by the amount of 56 Ni synthesized in the explosion, as well as the opacity of 
the ejecta (Shigeyama et al. 1992; Branch 1992; Khokhlov et al. 1993; Vacca & Leibundgut 1996; 
Hoflich et al. 1993, 1998; Dominguez et al. 2001). Changes in either of these, for example, due to 
changing progenitor metallicity with redshift, could affect the use of SNe la as standard candles. 
An evolutionary effect of 0.2 mag to z = 0.5 would nullify the SN la evidence that the Universe 
is accelerating, and measuring w to 10% requires that any effect be smaller than 0.04 mag. There 
are other routes to studying evolution with redshift: comparing SNe la in different host galaxy 
environments (Hamuy et al. 2000; Sullivan et al. 2003; Gallagher et al. 2005), and via detailed 
spectroscopic studies (Hook et al. 2005, Blondin et al. 2006, Bronder et al. 2006, in preparation). 
Neither approach has turned up any evidence of evolution. 

The rise-time has implications for SNe la explosion models. Varying the rise time from 20 to 
16 days at a fixed peak luminosity changes the implied amount of 56 Ni synthesized in the explosion 
by —10% (Contardo et al. 2000). Models of single white dwarf progenitor systems generally predict 
rise times in the range 13-19 days, while systems involving two white dwarfs allow for longer rise 
times because interaction with the disk of unaccreted material from the disrupted companion slows 
the diffusion of photons (Hoflich & Khokhlov 1996; Hoflich et al. 2002). 

The determination of the rise times has been the subject of some dispute. Historically, early- 
time, well-calibrated photometry of SNe la has been quite difficult to obtain. Some of the earliest 
studies of SNe la light-curve shapes examined the rise time (Pskovskii 1984). Riess et al. (1999a, 
hereafter R99) presented early observations of a set of nearby SNe la detected 10-18 days before 
maximum B luminosity. This data set was constructed primarily from unfiltered early detections, 
some from amateur observers. They transformed these observations to standard pass bands (in 
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particular, B) using models of early-time SN la spectra and colors. They then considered various 
late-time light-curve parameterizations, concluding that the rise time for a fiducial SN la was 
19.5 ± 0.2 days. Riess et al. (1999b) compared this number with the preliminary analysis of 
Goldhaber et al. (2001, hereafter G01) for distant SNe la and concluded that the rise times differed 
by 2.5 ± 0.4 days, a significance of ~ 6 a, and a clear signature of evolution. 

Aldering et al. (2000, hereafter AKN00) argued that this comparison was based on analyses 
that had ignored the significant correlations between the light-curve parameters, and that taking 
these properly into account increased the errors on the rise time to ±1.2 days for the distant sample. 
They concluded that the significance of the difference was closer to 1.5 a. 

Our purpose here is not to revisit this controversy; rather, we present a new, significantly 
more precise (~ 6 times) measurement of the high-redshift rise time and compare it with the value 
for nearby SNe. This has particular relevance because current SN projects place more stringent 
requirements on the standard-candle assumption. This paper presents measurements of the rise 
time from one such survey, the Supernova Legacy Survey (SNLS,Astier et al. (2006)). The design 
of this (and some other modern surveys) results in a dramatic increase in the amount of early-time 
photometry available when compared with previous generations of surveys, as detailed in §3. 

We first describe the basic problem (§2), then present the data on which our measurement 
is based (§3), followed by a description of our analysis procedures (§4). Finally, we compare our 
measurements against nearby SNe la, as well as between different subsets of our data (§5). 

2. RISE-TIME PARAMETERIZATION 

In order to measure the rise time of our sample we require a mechanism for combining data 
from different SNe la, correcting for the differences in light-curve shape and peak flux, and a model 
for the early-time flux as a function of time. Considerable effort and ingenuity have been devoted to 
developing techniques for parameterizing SNe la light curves near and after maximum light (Riess 
et al. 1996, G01, Guy et al. 2005). Here we follow the stretch method as described in, for example, 
G01. The flux as a function of time after the rise-time region is represented as 

/(t)=/o^(r), (1) 

where /o is the flux at maximum, and ip is some normalized flux template appropriate to the 
passband under consideration, r is the effective date defined by r = (t — tmax) /s (1 + z), where 
^max is the date of maximum flux in some arbitrary filter (usually B), s is the stretch, and z is the 
redshift. Conventionally, s = 1 is defined to represent an average SN la. Once s, /o, and imax are 
fited to the data, we can combine data from different SNe la by converting from observed epoch t 
to r and dividing by /o to normalize the flux values relative to each other. We must also apply a 
K-correction so that the resulting data are all expressed in the same rest- frame filter. Note that 
different observations from the same SN are correlated by this procedure, and that in addition 
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estimates of the light-curve fit parameters are generally quite strongly correlated. We must take 
both into account in our analysis. For our purposes we limit the fit to this model to the "core" 
light curve between — 10 < r < 35 days. The lower limit arises because we fit the rise-time model 
in this range, and we want to prevent the fit from suppressing unusual rise-time behavior. The 
upper limit arises because after this effective epoch the SNe la enter the so-called nebular phase, 
in which the stretch prescription no longer works (G01). 

In the rise-time region (r < —10), we follow earlier work (R99, AKN00, G01) in making use 
of a simple quadratic model 

f(t) = a(r + r r f (2) 

for r > — r r , and at earlier times. r r represents the rise time for the "fiducial" s = 1 template 
defined by ijj. Because the template can vary from analysis to analysis, it is critical that any 
comparison between the nearby and distant samples use the same one. For our purposes, a is 
purely a nuisance parameter. This approach implicitly assumes that stretch continues to work well 
at early times, as shown by G01 and R99. Our data set is not particularly well suited to investigating 
more complicated relations between the rise time and stretch, although we have slightly more to 
say about this in §5. Note that we do not require continuity between the fit to the rise-time region 
and the core light curve. 

Following R99, the quadratic form can be motivated by simple physical arguments: at early 
times SNe la should be hot enough that the standard BV pass-bands are in the Rayleigh-Jeans 
tail of the spectral energy distribution (SED), so / oc r 2 T oc v 2 (r + r r ) 2 T where v is the velocity 
and T the temperature. This, coupled with the fact that the measured colors and velocities do not 
change rapidly compared to the time since explosion in this regime, suggests a (r + r r ) 2 behavior. 
Ultimately, one of our goals is to test whether this model adequately describes the data. However, 
it has some clear limitations. In particular, one could always change the physical rise time in a 
fashion we could not detect by adding an initial period of non-t 2 behavior at very low, and hence 
undetectable, luminosities. 

Because much of the R99 data comes from unfiltered, very broad band observations, it can 
only be mapped to a single rest-frame filter independently. Hence, we are not yet in a position to 
measure the B and V rise times separately for the nearby sample, and thus have restricted this 
analysis to rest-frame B. 

3. DATA 

The SNLS 1 relies on data from the deep component of the 5 year CFHT Legacy Survey, using 
the square-degree MegaCam imager on the Canada- France-Hawaii Telescope (Boulande et al. 2003). 



1 http://www. cfht.hawaii.edu/SNLS/ 
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Repeat g'r'i'z' imaging is performed on four fields every three to four nights of dark time. SNLS 
is a "rolling" search, in which the same fields are searched repeatedly for variable objects. Each 
observation acts both as a potential discovery image for new SNe and as follow-up for candidates 
already discovered in the same field, allowing a considerable multiplex advantage. The primary 
goal of this program is to measure the average equation-of-state parameter (w) of the dark energy 
to 5% using ~ 500 — 700 high-redshift SNe la (see (Astier et al. 2006) for cosmological results from 
the first year of data) . 

The pioneering high-redshift SN programs (Riess et al. 1998; Perlmutter et al. 1999) typically 
had a gap between the initial reference and the follow-up image used to discover new SNe of 
approximately 1 month. This was designed to catch SNe around maximum luminosity and allow 
follow-up observations during dark time, but it resulted in large gaps in the early-time coverage. In 
fact, the majority of SNe la discovered in this fashion have no early-time detections. In a fraction 
of cases, however, SNe were discovered after maximum, or the redshift and light-curve shape were 
such that one or two early data points were obtained. This allowed G01 to measure the rise time 
at z ~ 0.5 to be 18.3 ± 1.2 days (statistical errors only). In contrast, in a rolling search the gap 
between images is typically a few observer-frame days, and it is usually possible to go back to 
previous images after the SN is discovered and measure the early-time flux. Furthermore, the gaps 
in the light curve due to the lack of observations during bright time are essentially uncorrelated 
with the light-curve phase. 

In this paper we consider only SNe la with spectroscopic type confirmation. SN types were 
determined using observations with the Gemini (Howell et al. 2005), ESO-Very Large Telescope 
(S. Basa et al. 2006, in preparation) and the W.M. Keck telescopes (R. S. Ellis et al. 2006, in 
preparation). We make use of data from the first 2 years of SNLS and only consider SNe at 
redshifts below 0.88. Above this redshift, our i' filter maps most closely to rest-frame U for a SN la 
at maximum, and so to measure the B rise time there we depend heavily on our z' observations. 
These have a considerably lower signal to noise ratio because the efficiency of our CCDs falls off at 
these wavelengths. Furthermore, the z' data suffers from considerable fringing because the CCDs 
are thinned, which results in additional scatter and calibration uncertainties. Finally, the z' cadence 
is not as frequent as in the other filters. The resulting light-curve fits would be dominated by the 
rest-frame [/-band data, while we seek to measure the rise time directly in the B band. While it 
would certainly be possible to include these data, even with this restriction the systematic error 
is already dominant in our rise-time determination. In future SNLS data sets the z' situation will 
improve, both because of changes to the fringe processing and because more z' observations are 
now obtained in each dark-time cycle. 

For our rise-time sample we require at least one rest-frame B observation in the rise-time region 
(—30 < r < —10 days) and one near peak luminosity (—6 < r < 6 days) to ensure that different 
light curves can be accurately normalized to each other using fa. We also require at least one rest- 
frame U, B, or V observation between 10 and 35 effective days in order to ensure that the stretch 
is well determined. These requirements eliminate roughly half of our SNe la; the combination of 
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our redshift range and the lunar cycle results in many SNe with data in the rise-time region not 
having observations near peak. However, the fraction of usable SNe for this purpose still far exceeds 
that possible with earlier data sets. The requirements also mildly select against the lower redshift 
portion of our SN sample, so the median redshift of our rise-time sample is slightly higher than our 
overall sample. This results in a sample of 73 spectroscopically confirmed SNe la in the desired 
redshift range. We note that these requirements are far more stringent than those necessary for a 
cosmological analysis. 

SNLS has two independent photometric pipelines, one based in Canada and the other in Prance; 
this paper makes use of photometry from the Canadian pipeline. The details of this pipeline, as 
well as the resulting photometry, will be presented elsewhere. Briefly, photometry is performed 
using a non-parametric PSF fit to subtracted images, for which the PSF is derived from a set of 
'PSF stars' shared across all images of the same field. The SN images are never re-sampled; rather, 
the reference SN-free images are transformed and PSF-matched to the data images. As in Astier 
et al. (2006), we rescale our photometric errors by ~ +25% to take into account the correlations in 
the reference image induced by these transformations. 

It is obviously critical that the rise times of the nearby and distant samples be computed using 
the same method. R99 provides nearby photometry for epochs r < — 10, so we supplement this 
with data from a number of other sources in order to allow a consistent determination of /o and 
s. The sources are provided in table 1. There are some potential complications with combining 
the data sets. In R99, the early photometry is given as magnitude relative to peak as a function of 
epoch relative to peak without specifying the peak magnitude or the date of maximum precisely. 
For most of the SNe this is not a serious problem, as the light curves around peak are sufficiently 
well sampled that both quantities can be determined relatively unambiguously. For some of the 
SNe, R99 used published early-time photometry, which allows the data to be tied together exactly. 
However, we experienced problems with four SNe, which we excluded from the sample. First, 
we were unable to locate any published photometry for SN1996by. For SN1996bv and SN1998ef 
the published photometry either has poor coverage or is of sufficiently low quality that we cannot 
tightly constrain the peak and date of maximum. In addition to these three, SN1996bo has only 
V-band observations available in the rise-time region, so we have excluded it from our B rise-time 
determination. We note that the excluded SNe would add relatively little weight to the fit, as each 
SN only has 1 or 2 observations in the rise-time region. 

To the R99 sample we add two more recent nearby SNe la with good early coverage: SN 2001el 
and SN 2002bo. The latter, in particular, adds considerably to the data sample. We then have 
eight nearby SNe la: SNe 1990N, 1994D, 1997bq, 1998aq, 1998bu, 1998dh, 2001el, and 2002bo. 
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Table 1. Nearby SNe la used in the rise-time measurement 



SN 


z 


Reference 


1990N 


0.0034 


1,2 


1994D 


0.0015 


1,3 


1997bq 


0.0094 


1,4 


1998aq 


0.0037 


1,5 


1998bu 


0.0030 


1,4 


1998dh 


0.0089 


1,4 


2001el 


0.0039 


6 


2002bo 


0.0042 


7,8 



References. - (1) R99, (2) 
Liraet al. (1998), (3) Richmond 
et al. (1995), (4) Jha et al. 
(2006), (5) Riess et al. (2005), 
(6) Krisciunas et al. (2003), (7) 
Benetti et al. (2004), (8) Krisci- 
unas et al. (2004) 
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4. ANALYSIS 

The primary complication we face in determining the rise time is in handling correlations. 
Particularly at high redshift, many SNe la have multiple observations in the rise-time region. In 
order to combine observations from different SNe la the data must be flux normalized, and the 
epochs must be shifted to reflect the date of maximum and the timescale divided by the stretch 
and 1 + z. The result is a set of data points with significant correlations in both dimensions 
(time and flux). Furthermore, the parameters determining the amount of correlation (the light- 
curve-fit parameters) are highly correlated themselves. Directly addressing this situation through 
the covariance matrix of the early-time data is not entirely trivial. We have chosen a different 
approach, which is to use a Monte Carlo technique to handle the correlations. We find that, for 
our data sample, their effects are comparable in size to the contribution from the measurement 
uncertainties of the data points. In other words, an analysis that ignored these correlations would 
underestimate the final error by approximately y/2. Worse, the resulting estimate for r r would be 
incorrect by several tenths of a day. 

In order to fit the rise-time data, we need to know the light-curve fit parameters (stretch, date 
of maximum, /o). We determine these using the light-curve template of Knop et al. (2003) and a 
time series of spectral templates descended from those of Nugent et al. (2002). In our fits the flux 
scales are allowed to float independently in all filters, with only the stretch and date of maximum 
held fixed between different filters (the difference between the date of maximum in different rest- 
frame filters is set by our template). Our fitting procedure predicts the SED of the SN on each 
observed epoch in physical units, which can then be converted to the B flux. We do not include 
data bluer than the rest-frame U band in our fits, since there are very few observations of nearby 
SNe la to constrain our model SEDs in this region. This means that the g' filter is not used at 
redshifts above 0.4, which constitutes the majority of our sample. The outputs of this procedure 
are best-fit values for each of the light-curve parameters as well as their correlations. 

With the light-curve parameters in hand, we combine the data as described in §2, applying 
the .fT-correction based on the model predicted SED. We then fit the rise-time model defined in 
equation 2 to the rest-frame B fluxes in the range —30 < r < — 10 days. We minimize the x 2 °f 
the model with respect to the data, fitting to the individual observations. For the SNLS data there 
are usually multiple observations of each SN in each filter on each epoch, so we can reject outliers 
due to, for example, unidentified cosmic rays by removing data points that disagree by more than 
3.5 a with other observations on the same night. This does not catch all outliers because for some 
nights there are only one or two calibratable observations in a given filter due to weather or other 
issues. We therefore also apply a 3.5 a outlier cut with respect to the model fit in an iterative 
fashion. This raises the possibility that interestingly discrepant SNe could be removed from the 
sample, but we find that relaxing or removing this cut has no effect on the final answer except to 
increase the \ 2 of the fit. Furthermore, a SN-by-SN investigation of data points that are removed 
by this cut gives no convincing examples of unusual rise-time behavior. As for R99 and AKN00, 
varying the upper limit of the data included in our parabolic fit (r = —10) has little effect on our 
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results as long as it is earlier than r < — 8 days. 

This, however, does not address the question of how to handle the induced correlations in the 
rise-time region. One approach would be to use the error properties of the individual photometry 
points to generate a random realization of our data sample, fit the core light curve for each SN 
as described above to determine the light-curve parameters, and then use these to perform the 
fit in the rise-time region to the combined data. After repeating this many times, the resulting 
distribution of r r could be used to find the rise time and its associated errors. A useful simplification 
arises because, for the purposes of the rise-time fit, the only part of the core light curve we are 
concerned with are the light-curve fit parameters, and so we can work in this space instead of 
directly with the photometry to handle the correlations. The contribution to the final error on r r 
is then split into two terms: that arising from the correlations between different points induced by 
the flux normalization and conversion to r, and that from the random photometric noise of each 
measurement. 

To handle the first term, we use the covariance matrix between the light-curve fit parameters 
for each SN to randomly generate a large number of realizations (typically 2000) of those parameters 
for each SN in our sample using standard techniques (James & Roos 1975). Changing the light- 
curve fit parameters changes the model SED for each data point, and so the K-corrections must be 
re-calculated and applied with each realization. For each of these sets we fit the rise-time region 
to determine specific values for r r and a, and then combine the results from all of the realizations 
to get the final values and their associated errors. The cost of this approach is that we cannot 
associate a simple \ 2 statistic with our overall fit. However, fits to individual realizations generally 
result in acceptable \ 2 values (for the main result presented below the \ 2 is I860 for 1409 degrees 
of freedom), especially considering that the correlations are not included in these numbers. This 
indicates that the quadratic rise-time model is a good representation of our data. The correlations 
affect both the fit value and errors for r r significantly. 

We can then include the term from the variances of the individual points by randomly re- 
sampling the rise time and a values using the error reported for each individual rise-time fit. The 
resulting errors on r r can be found by finding the ranges that contain the desired fraction of the total 
probability around the mean value. We use limits that are symmetric in probability space around 
the mean. The distributions have somewhat non-Gaussian tails, so the 2 a errors are generally not 
exactly twice the size of the the 1 a errors, etc. 

5. RESULTS 

For our sample of 73 high-redshift SNe la, we find the rise time to be r r = 19.101q 17 days (sta- 
tistical errors only). The histogram of rise-time values is shown in figure 1. This is approximately 
a factor of 6 more precise than the measurements of AKN00 or G01. The low-redshift SN sample 
gives T r = 19.58^o'i9- The nearby data in the rise-time region are shown in figure 2. We resist 
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comparing these until we have estimated the systematic errors below. The errors can be roughly- 
checked using a bootstrap analysis, which agrees with those values quoted above. Our estimates 
for the nuisance parameter 1000 x a are 6.I5IQ3} for the SNLS sample and 5.65l '27 f° r ^ ne near by 
one. a and r r are quite correlated, with a correlation coefficient of p = —0.6. 

In order to test the adequacy of the quadratic rise-time model, we have also performed fits in 
which the exponent is allowed to vary in the rise-time relation (n in / oc t n ). The best constraint 
comes from low redshift, at which we find n = 1.7±0.2, with the rise time reduced to r r = 18.80^Q32 
days. Not surprisingly, the errors on r r are considerably larger if n is not held fixed. The SNLS 
sample gives n = 2.0^Q3, with r r = 19.39^0 §2 days. The combined value is n = 1.8±0.2, essentially 
consistent with the assumed value of 2. In order to test that our conclusions are not too dependent 
on the value of n, we also re-fit the nearby and distant samples using a fixed value of n = 1.8, 
finding rise times of r r = lS^to'jJ and 18.491qJ 5 days, respectively. Fixing n = 1.8 shifts the rise 
time but does not appreciably affect the difference between the two samples. This conclusion also 
holds true when we compare subsets of the high-redshift sample, so we restrict the discussion to 
n = 2 subsequently. 

It is also quite interesting to compare the rise-time measurements for different subsamples of 
our data. The results of these fits are given in table 2. We do not consider these subsets of the nearby 
sample because it is too small. First, we consider splitting by redshift to search for evolution within 
our sample. The median redshift is 0.647. This is quite close to the transition between observer 
frame r' matching with rest-frame B and i! 1— > B, which takes place at z = 0.589. This test is 
therefore sensitive to two possible effects: evolution, and a calibration mismatch between r' and 
i'. Because we cannot test each of these independently, we split at z = 0.589, which results in 29 
SNe la in the intermediate- z sample and 44 in the high-z sample ((z) = 0.43 and 0.74, respectively). 
The photometric noise in the high-z portion of the sample is much larger than in the intermediate- z 
portion, as shown in figure 3. This figure also demonstrates that data from different SNe la can be 
combined quite accurately using the techniques described in §2. 

The data in the rise-time region are shown in figure 4, again split into the two groups. Using 
the same analysis, for z < 0.589 we measure r r = 19.01^Q'}g days and for z > 0.589 we measure 
r r = 19.67^Q49 days. The intermediate- z portion of the sample clearly dominates the fit to the full 
sample. These are statistically compatible (the difference is 1.2 a). 

As a test of whether the stretch model works at such early times, we split the sample by 
stretch. Unlike the nearby sample, one cannot measure the rise time precisely for most of the 
distant SNe la individually, so this test must be done on a sample-wide basis. Splitting around 
the mean stretch, the low-stretch sample (s < 0.99, (s) = 0.92) of 36 SNe gives r r = l^Oi^ 
days and the high-stretch sample (s > 0.99, (s) = 1.05, 37 SNe) gives r r = 19.09^Q2o days. These 
straddle the full sample value, and are statistically indistinguishable (0.3 a difference). If stretch 
did not work at early times, we generally would not expect these to agree. 

We can also split the sample by host galaxy star formation rate, following the methodology 




Fig. 1. — Rise-time (r r ) histogram for 73 SNLS SNe la, taking into account the correlations between 
data points as described in §2. This represents 2000 realizations of the light curve parameters. The 
mean value is 19.10 days. 
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Fig. 2. — Data in the rise-time region for the eight nearby SNe la used in this study. These have 
been combined using the techniques described in §2. 
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Fig. 3. — Shifted, normalized, and -fC-corrected SNLS data split by redshift and over-plotted. The 
rest-frame B band is shown. Observations in the same filter on the same night are combined for 
display purposes; the actual fit is performed on the individual data points. The blue triangles 
represent redshifts where rest-frame B matches the r' filter, and the red squares where i' is the best 
match. 
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Effective Day 



Fig. 4. — Data in the rise-time region for the 73 SNLS SNe la used in this study, combined using 
the techniques described in §2. They are split by redshift as in figure 3. Observations from the 
same filter on the same night have been combined for display purposes. 
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of Sullivan et al. (2006). Briefly, the broadband colors of the host galaxies are fitted by a galaxy 
spectral evolution code using the known redshift to determine the star formation rate per unit mass 
(sSFR), measured in yr~ x . The sample is then split between SNe in hosts with no star formation 
(passive, zero sSFR), those with moderate star formation rates (active, — 12 < log 10 sSFR < —9.5), 
and those with a large amount of star formation (vigorously star forming, log 10 sSFR > —9.5). We 
follow the above paper in limiting the application of this technique to z < 0.75, where it is most 
reliable for our data set. This results in a sample of 9 SNe la in passive hosts, 11 in active galaxies, 
and 35 in vigorously star forming galaxies. Clearly, this is an area where an increased sample size 
would be beneficial. Note that the comparison of these subsets is not independent of the stretch 
comparison: there is a relation between stretch and host galaxy morphology (Hamuy et al. 2000) 
and, in addition, star-formation rate (Gallagher et al. 2005) for nearby SNe la that has recently 
been confirmed at high redshift in the SNLS sample (Sullivan et al. 2006). In any case, for the 
passive sample we find a value of r r = 20.401{'5q days, for the active sample r r = 18.95^45 days, 
and for the vigorously star-forming sample r r = 19.07^Q'jg days. The rise time in passive hosts is 
mildly different than the others but only at the 1.2 a level. The data is shown in figure 5. 

The dominant systematic error for our measurement should arise from the ii'-corrections. 
Given the relative rarity of early-time photometry of nearby SNe la at early epochs, it is not 
surprising that there is very little spectroscopy available in the rise-time region. Hence, our ability 
to combine data from SNe la at different redshifts accurately is highly dependent on the theoretical 
models used to derive the SED at these epochs. At sufficiently early times SNe la should have 
SEDs dominated by thermal continuum. However, it is not observationally clear at what point this 
becomes a poor approximation. 

We have tried to quantify this by considering the effects of using a different, largely inde- 
pendently derived set of spectral templates, specifically those of Nobili et al. (2003). Redoing 
the above analysis, including the light-curve fits, shifts the rise time by 0.2 days. We take this 
as an estimate of the systematic uncertainty arising from if-corrections. AKN00 also considered 
the effects of changing the late-time light-curve behavior, arguing that this gave an upper limit 
for the systematic error of 2-3 days. However, the conditions considered were somewhat extreme, 
namely that all of the high-redshift sample have unusual late-time behavior when compared with 
the nearby sample. As an upper limit this is reasonable, but as a systematic estimate it is quite 
conservative. We expect that the systematic error involved in comparing different subsets of the 
SNLS sample against each other should be smaller, since uncertainties in the SED should affect the 
two samples in a similar fashion. In any case, since the subset comparisons reveal no effect even 
without systematic errors, we have not tried to estimate them. 

Ideally for the low-redshift sample .KT-corrections should not be a problem. The SNe la in the 
current nearby rise-time sample are all essentially at zero redshift, and so if the observations were 
on the standard photometric system there would essentially be no SED dependence of this process. 
Unfortunately this is not the case; a significant portion of the nearby data comes from unfiltered 
observations. We are not in a position to use the same technique to estimate the systematic error 
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Effective Day 

Fig. 5. — Shifted, normalized, and if-corrected SNLS data split by host galaxy specific star for- 
mation rate. In all other respects this is identical to figure 3. The blue triangles include SN from 
both active and vigorously star forming hosts as defined in §5 (log 10 sSFR > —12), while the red 
squares show those from from passive hosts. 



Table 2. Results of rise-time fits 



fit iVg N r r (days) a Notes 



SNLS 


73 


19.10l{u? 


0.42 
0.36 




nearby 


8 


19.581°;^ 


0.46 
0.38 




SNLS intermediate z 


29 


19 m +0 - 18 


0.37 
0.36 


z < 0.589, (z) = 0.43 


SNLS high z 


44 


19.67±g;S 


1.13 
0.99 


z > 0.589, (z) = 0.74 


SNLS low s 


36 


19.20l°;i 


0.67 
0.67 


s < 0.99, (s) = 0.92 


SNLS high s 


37 


io nq+ - 20 


0.39 
0.40 


s > 0.99, (s) = 1.05 


SNLS passive hosts 


9 


20.401^ 


1.73 
2.06 


log 10 sSFR < -12.0 


SNLS active hosts 


11 


i8.95tj}:j2 


0.74 
0.83 


-12 < log 10 sSFR < -9.5 


SNLS vigorous hosts 


35 


19.0718:11 


0.40 
0.37 


log 10 sSFR > -9.5 



a Both the 68.3% and 95.4% confidence limits are given (1 and 2 a, respec- 
tively) . 

Note. — Rise-time fits to various samples. sSFR is the star formation rate 
per unit mass in units of yr _1 . 



-18- 



on the nearby data. R99 discussed some of the uncertainties associated with transforming to the 
standard pass bands, but choose to include these effects as large statistical errors rather than as 
an overall systematic effect. Therefore, we simply have to trust that our derived statistical errors 
incorporate systematic uncertainties due to fT-corrections for the nearby sample. 

6. CONCLUSIONS 

We have measured the rise time from a sample of 73 high redshift (z = 0.15 — 0.9) spectroscopi- 
cally confirmed SNe la discovered and observed by the SNLS. This determination is roughly 6 times 
more precise than those previously available in this redshift range (AKN00, G01). Our measure- 
ment for this sample is r r = 19.10^q if (stat) ± 0.2 (syst) days. Using the same analysis technique 
on a sample of eight nearby SNe la [z < 0.1), we derive a value of r r = 19.58^g'^g days, where the 
quoted error incorporates both statistical and systematic errors. These differ at the 1.4 a level. In 
other words, using a considerably more precise comparison made possible by a substantially better 
data set, we find no compelling evidence for any difference between the rise times of nearby and 
distant SNe la. It is important to understand the limitations of this measurement in terms of its 
constraints on theoretical models. As was the case in R99, AKN00, and G01, the uncertainties 
presented above are the error in the mean stretch-corrected rise times of the two samples, not the 
scatter of rise times between individual SNe la. However, testing for differences between the two 
samples is still a very useful check against evolutionary effects that may be affecting cosmological 
analyses using SNe la. 

The above result suggests that we are currently limited by the systematic uncertainty associated 
with ET-corrections in performing this comparison. Therefore, significant advances will likely require 
better constraints on the early-time SEDs of SNe la. Alternatively, for a large enough sample, it 
may be possible to constrain the rise-time by only considering redshifts for which the observer and 
rest-frame filters match particularly well, minimizing the i^-correction errors. 

The quadratic rise-time model, motivated by simple physical arguments, provides a good fit 
to the data. Dropping this assumption, we find n = 1.8 ± 0.2 for / oc t n at early times. Fitting 
for this extra parameter substantially weakens our constraints on r r , but does not indicate any 
discrepancies. Redoing the fits with n fixed at this value also does not significantly affect any of 
our results, simply decreasing all of the measurements of r r by ~ 0.6 days. 

We have also split the SNLS sample into subsets and searched for differences between them. 
The comparisons we consider are splitting the sample by redshift at z = 0.589 to test for evolution 
within our sample and for calibration systematics, splitting by stretch, and splitting by host galaxy 
star formation rate. In all of these cases the subsamples give compatible rise times. The fact that 
the low- and high-stretch samples agree confirms the claim of G01 that the stretch parameterization 
works both at early and late times, at least up until around day 35. Unlike the rise-time fit to 
the full sample, for many of these subsets we are limited by statistics, so as the survey continues 
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more sensitive comparisons will be possible. Once a full cosmological analysis on all these SN is 
available, it will also be interesting to check whether there is any correlation between the residual 
from the best fitting cosmology and rise time. 
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